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Abstract 



We discuss the effects of the s-channel unitarization on the x and dependence of 
structure functions. The unitarization is implemented at the level of photoabsorption 
cross sections by resorting to the light-cone wave functions of virtual photons and to 
the diagonalization property of the scattering matrix in a basis of Fock states of the 
photon with fixed transverse size. Triple pomeron effects are also explicitly taken into 
account. We find large unitarity corrections to the structure functions at x < 10~^. 
The results are in very good agreement with the existing NMC and the preliminary 
HERA data. 
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1 Introduction. 

The conventional QCD evolution predicts a rapid rise of the density of partons in the 
Regge limit of deep inelastic scattering, l/x 1 and large (for a review see 
The gluon distribution is predicted to increase as 

xg{x, Q2) oc exp \c^\og [l/a^lQ^)] log (1/x)} . (1) 

and the same behaviour is expected for the sea density. This leads to violation of the 
s-channel unitarity and of the Froissart bound at large s = Q'^/x. 

L. V. Gribov, Levin and Ryskin (GLR) proposed a modified nonlinear evolution 
equation to account for the recombination of partons which suppresses the increase of 
their densities in the region of very small x 0. 

In this paper we suggest to look at the unitarity problem from a different point of 
view, working at the level of cross sections. The starting point is that the photoab- 
sorption cross section can be derived as an expectation value of the interaction cross 
sections of the multiparton Fock states of the virtual photon over the wave functions 
of such states 0, The s-channel unitarity must be imposed on the partial- wave 
amplitudes of interaction of these Fock states of the photon, and this unitarization of 
the partial-wave amplitudes can be reformulated in the parton model language as the 
fusion of partons. As far as the x and dependence of the total photoabsorption cross 
section is concerned, the approach Q is equivalent to the conventional Weizsacker- 
Williams reinterpretation of deep inelastic scattering in terms of the (tree-level) 
Compton scattering off the radiatively generated sea quarks, but has few technical ad- 
vantages: i) an exact diagonalization of the scattering matrix in the representation of 
the Fock states of the photon, ii) an easy identification of the partial-wave amplitudes 
which are subject to the s-channel unitarization. 

Understanding the unitarization effects becomes crucial in view of the forthcoming 
data from the HERA experiments, which will go down to x ~ 10~^. In this paper we 
shall offer a semiquantitative analysis of the unitarity effects at large 1/x, including 
the estimate of the triple pomeron (fan diagram) contributions, the driving term of 
which arises from the qqg Fock component of the photon. We extend to small x, 
including the unitarity effects, predictions for the unitarized structure function of the 
proton F2{x, Q^), based on the dynamical generation of the gluon and sea distributions 
starting with the three constituent quarks |^ 

2 The color dipole cross section as a subject of uni- 
tarization 

At 1/x ^ 1 photoabsorption can be treated as the scattering on the nucleon of the 
qq pairs the virtual photon transforms into at large distances Az ~ 1/m^x ^ 
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upstream the target nucleon (Fig. 1). In this hmit the transverse size p of the qq pair 
and the Sudakov variable z , i.e. the fraction of the photon's hght-cone momentum 
carried by one of the quarks of the pair (0 < 2; < 1), are frozen in the scattering process. 
The scattering matrix becomes diagonal in the (p, z) representation, which allows one 
to write the transverse (T) and longitudinal (L) virtual photoabsorption cross sections 
as the quantum mechanical expectation value 

(TtAQ^) = I'd^ I d'P |*t,l(^,p)|Vo(p) . (2) 

Here cro(p) is the interaction cross section for the qq color dipole of size p, given by 

16 f d^kV{k)[l- exp{ik- p)] ^ 



Mp) = y«.(p) J ^^-^^ asik ) . (3) 

In k is the transverse momentum of the exchanged gluons (Fig. la), /^g ~ l/i?c is 
the effective mass of gluons introduced so that color forces do not propagate beyond the 
confinement radius Re- The gluon-gluon-nucleon vertex function V{k) = 1 — J^chi^k"^) , 
where Tchi^f) is the charge form factor of the proton. 

The wave functions of the qq fluctuations of the photon were derived in ||^ and read 

\^t{z.p)? = ^E4{[-' + (1 - ^fVK^i^Pf + mjMepn , (4) 

\^l{z,p)\' = ^'(1 - ^fM^pr > (5) 

where Ky{x) are modified Bessel functions, = z{l — z)Q'^ + m'j and is the quark 
mass. Since either V^(k) and/or the factor [1 — exp(ik- p)] vanish as |k| — >■ 0, the dipole 
cross section cro(p) is free from infrared divergences. At small p 

o-o(p) oc asip) log [1/ asip)] , (6) 

and vanishes as p ^ 0. 

Eq. d^) gives the driving, constant term of the QCD pomeron contribution to the 
dipole cross section. The x dependence of the dipole cross section ao{x,p) comes from 
the t-channel iteration of the gluon-exchange between the two gluons in the diagram of 
Fig. lb. These higher order contributions can be understood as the QCD evolution of 
the density of gluons in the nucleon, or interactions of the higher qqgi...gn Fock states 
of the photon [0. The convenient way to incorporate this x-dependence in is as follows: 

The driving term G{x, Q"^) of the distribution of soft, a; <^ 1, gluons with < k^ < 
Q^, generated perturbatively from the nucleon, is given by 0, ^ 

G{x, Q') = -- asik') V{k) . (7) 



TT X Jo (k^ + p|,) 



3 



Comparing eq. @ and eq. (0) we can rewrite the dipole cross section as 



Mp) = — as(p)y^(l-e-^)-4^^. (8) 

The generahzed QCD ladder diagrams of Fig. lb are summed by replacing the driv- 
ing term of the gluon flux G{x, Q"^) with the full GLDAP-evolved gluon distribution 
giGLDAP^x^Q"^) of our previous work |]^, i.e. 

. , 47r2 /"dk^ ■^,Q^dxg{GLDAP,x,k'^) 

This GLDAP dipole cross section (y) must then be used in calculation of the pho- 
toabsorption cross section and, eventually, of the GLDAP structure function using 
the definition F2(a;,Q^) = Q^cr^.7v(a;, Q^)/47r^aem- Since o"(x,p) ~ (?as{p)g{x-,Q'^\ 
in view of Eq.(||) at small x it violates both the Froissart bound and the s-channel 
unitarity. 

Now, let us identify the corresponding partial waves using the conventional impact 
parameter representation for the scattering amplitude 

/o(q) = 2z 1 d^h ro(b) exp(-?b • q) . (10) 

For the predominantly imaginary scattering amplitude /o(q) = '^o"o(x, p) exp(— i?oq^/2) 
the profile function (partial- wave amplitude) ro(b) of the GLDAP dipole cross section 
(9) equals 

and violates the s-channel unitarity constraint r(b) < 1 when aQ[x^p) becomes large. 
(By geometrical considerations, the diffraction slope can be related to the conven- 
tional hadronic slope for vrp scattering through Bq = -B,rp/2 + l^)- 

We encounter, in fact, the as yet unsolved problem of unitarization of the rising 
strong interaction cross section. Two different schemes are commonly adopted for 
construction of the unitarized partial-wave amplitude r(b) : the eikonal unitarization 



m 



r(b) = 1 - exp [-ro(b)] , (12) 
and the K, matrix unitarization [0 

1 + ro(b) 

In both cases as the bare profile function ro(b) rises, the unitarized one r(b) tends to 
the black disk limit r(b) 1. 
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The unitarized total cross sections cr{x,p) =2/ cPhT{h) reads 

eikonal : cr(x, p) = AnB [\og{ri{x, p)) + Eiirj) + 7] ; (14) 

/C matrix : cr(a;, p) = AnB log(l + ri{x, p)) , (15) 

where 7 is the Euler-Mascheroni constant, Ei is the integral exponential function. The 
quantity ri{x, p) = cro{x, p)/47ri?o shown in Fig. 2, controls the effect of the unitarization: 
a{x,p) ~ ao{x,p) at ri{x,p) <^ 1, and at ri{x,p) ^ 1 the unitarization suppresses the 
cross section, a{x,p) <^ cro(a;,p). This plot shows that the unitarization effects are 
important only at large p. The difference between the two cross sections (|T^) and 
([T5|) can be taken as a measure of the theoretical uncertainty on the unitarization 
procedure. For convenience we shall adopt hereafter the /C-matrix procedure, the 
eikonal unitarization gives similar, if somewhat weaker unitarizations effects. 



3 Unitarity and the triple pomeron contribution 

To the leading order in the s-channel unitarization r(b) = ro(b) — xro(b)^/2, where 
X = 1(2) for the eikonal (K- matrix) unitarization, and the unitarized photoabsorption 
cross section reads 

(TrN{x, Q2) = 2 1 d^h r(b) = (ao) - x(/ d^hV,{hf) , (16) 

where we have denoted by (•) the average over the photon wave function. The first 
term in eq. (|IBp gives the GLDAP cross section and the conventional GLDAP structure 
function, the second term is the unitarity (shadowing) correction, which in principle 
is a nonlinear functional of the GLDAP cross section. The unitarity correction can be 
related to the cross section of the forward diffraction dissociation of the virtual photons 
7* + p — > X + p, where X = qq, as follows : 



danil* qq) 





(17) 



dt 



dMMt 



i=0 

Here Bd stands for the diffraction slope of the diffraction dissociation. The driving 
term of diffraction dissociation is excitation of the qq component of the photon, which 
has the oc 1/{Q^ + M'^f mass spectrum and gives the x-independent shadowing term. 
Diffraction excitation of the higher qqgi---gn Fock states of the photon gives rise to the 
triple pomeron component of the mass spectrum at 3> Q^: 



daoil* qqg...) 



dtdM2 



= G31P ^ , (19) 

i=0 



5 



and to the logarithmically, oc log(s) oc log{l/x) , rising shadowing component. 

We shall present a brief discussion of the effects of the qqg Fock state, which is the 
driving term of the triple-pomeron mass spectrum. We apply to the qqg system the 
same wave function formalism already used for the qq component: 

(20) 

Here we have introduced the three-particle wave function $ and the three-particle 
interaction cross section o"o(r, R, p), where r, R and p are the gq, qg and qq separations, 
respectively, and Zg is the fraction of the photon's light-cone momentum carried by the 
gluon. It is possible to show that the leading contribution to the triple pomeron 
component comes from the ordering l/Q^ -C -C ~ i?^. In this case the gg-pair 
of the qqg Fock state can be treated as the pointlike color-octet charge, so that the 
three-particle cross section ao{r,R,p) reduces to 

cro(r,i?,p)^^ao(r), (21) 

where 9/4 is the familar ratio of the octet and the triplet strong couplings. The three- 
particle wave function takes on the factorized form 

I^P = TT^-«.(p)l^(^,p)p4-^(/^Gr), (22) 

which holds when the gluon of the qqg Fock state participates in the interaction, so 
that only the diagram of Fig. Ic (and not the one of Fig. Id) should be taken into 
account and the corresponding three-particle cross section is reduced to half of (|21|), 
we refer to |p for a derivation and more details. Here \E'(2;,p) is the wave function of 
the qq state and J-'{x) = {x'^K[{x)Y satisfies ^(0) = 1 and J-'{x) ~ exp(— 2a;) at large 

X. 

Since at ^ one has dM'^/M'^ = dzg/zg, the 1/zg dependence of the wave 
function (^2]) , which has its origin in the spin 1 of gluons, leads to the mass spectrum 
of the form (19) and to the oc log(l/a;) dependence of the shadowing term. Using the 
wave function (22) and the cross section (21) it can be shown that the triple-pomeron 
component of the diffraction dissociation cross section (20) satisfies an approximate 
factorization property: 



dtdM2 



(Tj^nAstp— (23) 

=0 



where A^jp only weakly depends on x, and on the flavour of quarks in the qq 
state. This approximate factorization allows one to write down the unitarized structure 
function in the form 

F^ix, Q') = F^iGLDAP, x, Q') (l - ASh{qq) - log(l + ^)) . (24) 
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Here we singled out the shadowing correction from excitation of the qq Fock state, 
which has a strong flavour dependence, ASh{qq) oc l/m'j, and vanishes oc l/Q^ for the 
longitudinaly polarized photons [P, |, |12|. In Eq. (24) B^jp stands for the diffraction 
slope of the diffraction dissociation in the triple-pomeron region, and xq ~ 0.1 corre- 
sponds to the usual definition of the diffraction excitation as interactions in which the 
target nucleon receives little recoil and emerges in the final state separated from the 
virtual photon's debris by the rapidity gap At] = log(l/xo) > 2. Diffraction excitation 
of the qq Fock state dominates the mass spectrum at ~ and is the counterpart of 
excitation of resonances in the hadronic scattering, which has the diffraction slope 
close to that of the elastic scattering, so that here we take Bd = 10 {GeV/c)~'^ ~ B.„n, 
whereas in the triple-pomeron region Bd = B^jp ?a i?^7v/2 (for the review see |1T3|). 
More complete treatment of the shadowing, in which one describes the diffraction dis- 
sociation in terms of the structure function of the pomeron, does not change much the 
above simple estimate of the triple-pomeron effect. 



4 Discussion of the results 

In our approach to DIS at small-x the fundamental quantity is the dipole cross section 
(To(x, p). We assume that at large p ~ Rn it is practically fiat as a function of x, 
which gives a smooth connection with the slow rise of the hadronic cross sections 
I^ . This corresponds also to our treatment of the glue in nucleons at < Ql as 



a radiative effect, starting from three valence quarks at the spectroscopic scale, as 
explained in detail in our previous papers ^ |^. Beyond Ql we calculate g{x,Q'^) 
using the conventional GLDAP evolution. At small we use the freezing coupling 
constant a<i(Q^ < Q/) = C(s{Q}), where Q'j is a scale of the same order in magnitude as 
Ql- Also, at < Q2 

we neglect the splitting of gluons and the momentum flow from 
glue to sea. The absolute normalization of the sea structure function is determined 
by the absolute normalization of the dipole cross section cro(p), which is fixed by the 
confinement radius Rc = l/pc and the frozen strong coupling asiQ})- These two 
parameters pc and Q / are strongly correlated when we ask for simultaneous description 
of the ttN total cross section. This procedure corresponds to the minimal infrared 
regularization and has been already applied with success to various structure function 
calculations [|], ^ ||, |12|- The choice of parameters used in the present calculations is: 
Ql = 1.0 GeVyc^, Q) = O.SGe^Vc^, pc = 300 MeV/c\ 

With the above set of parameters, in the region of small we find A^jp ~ 
0.15 GeV""^ , to be compared with the experimental result for the diffraction disso- 
ciation of real photons A^jp = 0.16 GeV~^ |T^. As stated above, our calculations show 
that A^jp changes little, < 10%, with and on the flavour of quarks. 

We predict the GLDAP and unitarized structure functions with absolute normal- 
ization. The predictions for the unitarized structure function F2{x,Q'^) at different x 



and Q are shown in Fig. 3 and Fig. 4. compared to the existing NMC data |]I5| at 
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X = 0.025 and to the preliminary HERA data fl^, [H at a; = 1.5 ■ 10'^, A - 10'^. The 
unitarity correction is numerically large and makes the rise of the structure function 
with 1/x much less steep. Considering the accuracy of the NMC data, the shadowing 
correction is rather large already at x ~ 10~^, where the effect of shadowing of the 
qq Fock state of the photon and the triple-pomeron contribution are of comparable 
magnitude. At smaller x the triple-pomeron component starts taking over, but the qq 
component also rises slowly as it is proportional to pY). We find very good agree- 
ment between our predictions for the unitarized structure function and the experiment. 
Notice, that the triple pomeron component tends to dominate as x decreases. 

In Fig. 5 we compare our prediction for F2 in the very small x region with some 



parametrizations [|T^, based on fits to the previously published experimental data 
on structure functions. At very small x ~ 10~^ our structure functions follow an 
approximate law F2(x, Q^) oc (l/x)^^'^ \ where the exponent A((5^) ~ 0.21 at = 
AGeVyc^ and A{Q^) ^ 0.311 at = IhGeV^/c^. 

In Fig. 6 we present the effects of the unitarization on the charm component of 
the structure function. We take rric = 1.7 GeV/c^. Here the shadowing of the qq Fock 
state is negligible, and the unitarity correction is completely dominated by the triple- 
pomeron component. The situation with the ratio R = a ^ / (Jt is similar: the ASh{qq) 
component is small, and the unitarization of the longitudinal structure function is 
dominated by the triple pomeron component. We find that R is quite structureless and 
depends very weakly on at very small x. At x < 10~^ we predict R ~ 0.20 — 0.24. 

After this work was completed, there appeared an estimate of the unitarity effects 
by Askew et al. [2C], who use a very different formalism. In |20| the parameter which 
controls the shadowing is the radius Rg of the region of the proton in which the gluons 
are concentrated. With Rg = 2GeV~^, which is significantly smaller than the proton's 
size. Askew et al. find unitarity corrections close to ours at x < 10"'^, and have no 
shadowing at x ~ 10~^. 



5 Conclusions 

We have presented a simple approach to the unitarization of rising structure functions 
at small x. The shadowing term in the unitarized structure function is dominated by 
the triple-pomeron term, which is approximately flavour- and Q^-independent. We 
find large shadowing corrections already at relatively large x ~ 10-2. We find that the 
dynamical generation of the glue and sea from the valence quarks, which predicts the 
small-x structure functions with absolute normalization, is in good agreement with the 



preliminary data from HERA WTl [IB . 
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Figure captions: 



Fig.l - a) Driving term of the QCD pomeron, b) Generalized ladder diagrams for the 
QCD pomeron, c) Interaction of the qqg Fock state of the photon which gives the 
rising contribution to the total cross section, d) Interaction with the color-octet 
qq component of the qqg state which does not contribute to the cross section 
growth. 

Fig.2 - The unitarization parameter ri{x,p). 

Fig. 3 - Predictions for the GLDAP (dashed curve) and unitarized (solid curve) structure 



functions at small x compared to the NMC [jl5[ data at a; = 0.025 and the 
preliminary HI [|l^] and ZEUS [jl6l data at small x. The dot-dashed curve shows 
the effect of shadowing of the qq Fock state of the photon. 

Fig.4 - Our prediction for the x dependence of the unitarized structure function at 
= 15 (GeV/cf compared to the NMC 



HI ig and ZEUS [|T|] data. 



Fig.5 



Fig.6 



- Comparison of our prediction for the small-x behaviour of F2{x, Q"^) (solid curve) 
with some of the recent MRS and DGJLP [|T^ parametrizations based on 
fits to the published data. 

- Our prediction for the GLDAP (dashed curve) and the unitarized (solid curve) 
charm distribution at = 15 {GeV/c)"^. The effect of shadowing of the cc Fock 
state of the photon is shown by the dot-dashed curve. 
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